A proposed search for a fourth neutrino with a PBq anti-neutrino source 
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Several observed anomalies in neutrino oscillation data can be explained by a hypothetical fourth 
neutrino separated from the three standard neutrinos by a squared mass difference of a few eV'^ . We 
show that this hypothesis can be tested with a PBq (ten kilocurie scale) ^"^^Ce or ^"^Ru antineutrino 
/3-source deployed at the center of a large low background liquid scintillator detector. In particular, 
the compact size of such a source could yield an energy-dependent oscillating pattern in event spatial 
distribution that would unambiguously determine neutrino mass differences and mixing angles. 



Most results from neutrino experiments over the 
last twenty years can be quite accurately described 
by a model of oscillations between three v flavors (z/e, 
Vfj,, Vr) that are mixtures of three massive neutrinos 
{vi, V2., vz) separated by squared mass differences of 
Aml^ = 8 • 10"^ eV^ and Am^^ = 2.4 • IQ-^ eV^ Q. In 
the past the LSND experiment suggested the existence 
of a fourth massive neutrino with a mass of ~1 eV^ [5j- 
This evidence has been confirmed by the MiniBoone 
experiment in the antineutrino sector [6[. Recently the 
hypothetical existence of a fourth v has been revived 
by a new calculation (2| of the rate of Ve production 
by nuclear reactors that yields a v flux about 3% 
higher than previously predicted. This calculation then 
implies [3[ that the measured event rates for all reactor 
Ue experiments within 100 meters of the reactor are 
about 6% too low. The deficit can also be explained by a 
hypothetical fourth massive v separated from the three 
others by |Ato^(,^| > 0.1 eV^. This mixing can explain 
a similar deficit in the rate of u interactions in Gallium 
solar-i/ detectors when exposed to artificial ^^Cr and 
•^^Ar MCi sources [J|. Combination of [3| and [4| deficits 
is significant at the 99.8% C.L., though no conclusive 
model can explain all data [7|. In this letter we propose 
an unambiguous search for this fourth neutrino by using 
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a 1.85 PBq (50 kCi) antineutrino source deployed at the 
center of a kilo-ton scale detector such as Borexino [9(), 
KamLAND [l3|, or SNO-I- [11]. Antineutrino detection 



FIG. 1. Simplified decay scheme of the pair ^''*Ce-^**Pr. 



will be made via the inverse beta-decay (IBD) reaction 
Ve + V ^ e++n. The delayed coincidence between de- 
tection of the positron and the neutron capture gamma 
rays will allow for a nearly background free experiment. 
The small size (~10 g) of the source compared to the 
size a nuclear reactor core may allow the observation of 
the characteristic !/-oscillation pattern of event positions. 



Large liquid scintillator (LS) detectors, called LLSD 
hereafter, share key features well suited to search for an 
eV-scale fourth v {p) state mixing with v^ {ve)- The 
active mass is composed of about thousand tons of ultra- 
pure LS contained in a nylon or acrylic vessel. The scin- 
tillation light is detected via thousands of photomulti- 
pliers uniformly distributed on a stainless steel spher- 
ical vessel. In Borexino and KamLAND the target is 
surrounded by mineral oil or scintillator contained in a 
stainless steel vessel. This buffer is enclosed in a wa- 
ter tank instrumented by photomultipliers detecting the 
Cherenkov light radiated by cosmic muons. In the follow- 
ing we study the deployment of a z/ source of energy spec- 
trum S(E^), mean lifetime r, and initial activity Ao, en- 
capsulated inside a thick tungsten (W) and copper (Cu) 
shielding sphere, at the center of a LLSD. We consider a 
running time te with a fully efficient detector. The theo- 
retical expected number of interactions at a radius R and 
energy E^ can be written: 



f^^l^ = Ao-n-a{E.yS{E.yV{R,E,) f\~''^dt, 
anatji, Jq 

(1) 

where n is the density of free protons in the target for 

inverse beta decay, a is the cross section. V{R, E^) is the 



2-1/ oscillation survival probability, defined as: 



V{R,E^) = l-sin-'(26'„eto)-sin^ 1.27 



Aml,JeV']R[m] 
£;^[MeV] 



(2) 
where Am^^,^ and ^now are the new oscillation param- 
eters relating Ve to the fourth ly. In our simulations 
we assume a 15 cm vertex resolution and a 5% energy 
resolution. In the no-oscillation scenario we expect a 
constant v rate in concentric shells of equal thickness 
(see Eq. [T]). At 2 MeV the oscillation length is 2.5 m for 
Am^(,„= 2 eV^, proportional to l/Am^g^ (see Eq. ^. 
A definitive test of the reactor antineutrino anomaly, in- 
dependent of the knowledge of the source activity, would 
be the observation of the oscillation pattern as a func- 
tion of the v interaction radius and possibly the v energy. 

Intense man-made v sources were used for the cal- 
ibration of solar-i^ experiments. In the nineties, ^^Cr 
(~750 keV, Ao -MCi) and ^^Ar (814 keV, ^=0.4 MCi) 
were used as a check of the radiochemical experiments 
Gallex and Sage [l6[. There are two options for de- 
ploying ly sources in LS: monochromatic z/g emitters, 
like ^^Cr or ^^^Ar, or D^ emitters with a continuous 
/3-spectrum. In the first case, the signature is provided 
by Ve elastic scattering off electrons in the LS molecules. 
This signature can be mimicked by Compton scattering 
induced by radioactive and cosmogenic background, or 
by Solar- ^' interactions. The constraints of an exper- 
iment with Vf, impose the use of a very high activity 
source (5-10 MCi) outside of the detector target. In the 
second option, Vf, are detected via inverse beta decay. 
Its signature, provided by the e~^-ii delayed coincidence, 
offers an efficient rejection of the mentioned background. 
For this reason, we focus our studies on Dg sources. 

A suitable De source must have Qp >1.8 MeV (the IBD 
threshold) and a lifetime that is long enough (>1 month) 
to allow for production and transportation to the de- 
tector. For individual nuclei, these two requirements 
are contradictory so we expect candidate sources to 
involve a long-lived low-Q nucleus that decays to a 
short-lived high-Q nucleus. We identified four such 
pairs "4ce-"4Pr (Q;3(Pr)=2.996 MeV), lo^j^^.ioej^j^ 
(Q^(Rh)=3.54 MeV), ^o^v-^Oy (Q^(Y)=2.28 MeV), and 
42Ar-42K (Qb (K)=3.52 MeV), some of them also re- 
ported in [l7|. The first three are common fission prod- 
ucts from nuclear reactors that can be extracted from 
spent fuel rods. While not minimizing the difficulty of 
doing this, the nuclear industry does have the technol- 
ogy to produce sources of the appropriate intensity, at 
the ppm purity level. In fact, 10 kCi ^'^Sr sources have 
been produced and used industrially for heat generation. 
Delays obtaining authorizations for transportation and 
deployment of the source into an underground labora- 
tory should be addressed at the start of the project. 
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FIG. 2. Advantage of i^e sources providing both R and Evis 
oscillation patterns. IBD rate for a 50 kCi ^*^Ce source de- 
ployed at a center of a LLSD, in 10 cm radius bins and 100 keV 
bins of visible energy, E„iB=Ee-|-2nie. In one year, 38,000 Ue 



and sin^(2Sn 



0.1. 



For this paper, we concentrate on the ^'*^Ce source 
(Fig. [T]) because its Qp is greater than that of ^"Sr and 
because it is easier to extract chemically than ^"^Ru. We 
note also that it has a very low production rate of high- 
energy 7 rays (> IMeV) from which the De detector must 
be shielded to limit background events. Finally cerium 
is present in fission products of uranium and plutonium 
at the level of a few percent. 



We now focus on the unique oscillation signature induced 
by an eV-scale sterile v at the center of a LLSD. For 
"4Ce-"4Pr, 1.85 PBq (50 kCi) source is needed to reach 
40,000 interactions in one year in a LLSD, between 1.5 
and 6 m away from the source (n^f =5.3 • 10^^ m^"^). This 
is realized with 14 g of ^^''Ce, whereas the total mass 
of all cerium isotopes is ~1.5 kg, for an extraction from 
selected fission products. The compactness of the source, 
<4 cm, is small enough to be considered as a point-like 
source for Arri^^^ eV^ oscillation searches. This source 
initially releases ~300 W, and it could be cooled either 
by convective exchanges with the LS, or via conduction 
though an ultrapure copper cold finger connecting the 
massive passive shield to a low temperature bath. /3^- 
decay induced Ve are detected through IBD. The cross 
section is a{Ee) = 0.956 10"'''^ x PeE^ cm^, where Pe and 
Ee are the momentum and energy (MeV) of the detected 
e^ , neglecting recoil, weak magnetism, and radiative 
corrections (%-level correction). The e"*" promptly 
deposits its kinetic energy in the LS and annihilates 
emitting two 511 keV 7-rays, yielding a prompt event, 
with a visible energy of Ee= Ei^-(?7i„-TOp) MeV; the 



Source |F.Y. "''"U/^^'^Pu| t^/j 1 1°' /3~ (keV)|2"'' /3~ (keV) | J.,>iAf,v 



'7>2MeV 



W/kCi|kCi/410^ int./y 



144cg_144pj. 


5.2%/3.7% 


285 d 


318 (76%) 
184 (20%) 
238 (4%) 


2996 (99%) 
810 (1%) 


1380 (0.007%) 
1489 (0.28%) 


2185 (0.7%) 


7.47 


43.7 


"'«Ru-i"6j^h 


0.5%/4.3% 


373 d 


39.4 (100%) 


3540 (78%) 
3050 (8%) 
2410 (10%) 
2000 (2%) 


1050 (1.6%) 
1128-1194 (0.47 %) 
1496-1562 (0.19 %) 
1766-1988 (0.09 %) 


2112 (0.04%) 
2366 (0.03%) 
3400 (0.016%) 


8.40 


23.0 



TABLE I. Features of ^""Ce-^^Pr and ^°''Ru-^°'^Rh pairs, extracted from spent nuclear fuel. F.Y. are the fission yields of "''Ce 
and ^"^Ru, ti/2, their half-lives. /3-end-points are given for 1°* and 2""* nucleus of each pair. The J-y's are the branching ratio of 
gammas 7 rays per beta-decay above 1 and 2 MeV. The two last columns are the heat produced/kCi and the activity required 
to get 40,000 events/year. 



emitted keV neutron is captured on a free proton with a 
mean time of a few hundred microseconds, fohowed by 
the emission of a 2.2 MeV deexcitation 7-ray providing 
a delayed coincidence event. The expected oscillation 
signal for Am'^^^ — 2 eV^ and sin^(20now) = 0.1, is 
shown on Fig. [2] LLSD are thus well suited to search for 
an eV-scale fourth ly state. Note that a study of signals 
of a ^"^Sr MCi source external of a LLSD was done in [19| . 

The space-time coincidence signature of IBD events en- 
sure an almost background-free detection. Backgrounds 
are of two types, those induced by the environment or 
detector, and those due to the source and its shielding. 

The main concern is accidental coincidences between 
a prompt (E>0.9 MeV) and a delayed energy deposi- 
tions (E>2.0 MeV) occurring within a time window taken 
as three neutron capture lifetimes on hydrogen (equiva- 
lent to about 772 /isec), and within a volume of 10 m^ 
(both positions are reconstructed, this last cut leading to 
a background rejection of a factor 100). The main source 
of detector backgrounds originates from accidental coin- 
cidences, fast neutrons, and the long-lived muon induced 
isotopes ^Li/^He and scales with R^ when using concen- 
tric i?-bins. These components have been measured in- 
situ for the Borexino geo-i^ searches [10[, at 0.14±0.02 
counts/day/ 100 tons. Being conservative we increase it 
to 10 counts/day/ 100 tons in our simulation. 

Geologic Pe arising from the decay of radioactive iso- 
topes of Uranium/Thorium in the Earth have been mea- 
sured at a rate of a few events/ (100 ton.year) in Kam- 
LAND [131 and Borexino [10]. Reactor D^ emitted 
by the /3-decays of the fission products in the nuclear 
cores have been measured in KamLAND at a rate of 
~10 events/(100 ton.year) in the energy range of inter- 



est [12|. We use a rate of 20 events/(100 ton.year), which 
is negligible with respect to the 9^ rate from a kCi source. 
The most dangerous source background originates 
from the energetic 7 produced by the decay through 
excited states of ^'*''Pr (Table |T|. We approximate 7 
ray attenuation in a shield of 33 cm of W and 2 cm of 
Cu with an exponential attenuation law accounting for 
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FIG. 3. Signal rate of a 50 kCi ^''''Ce deployed for 1 year at a 
center of a LLSD (black), compared to the sum of all identified 
backgrounds rate (blue), as a function of the detector radius 
in 10 cm concentric bins. A shield made of 33 cm of W and 
2 cm of Cu attenuates the backgrounds, dominated primarily 
by ^^*Ce 7 lines, then by external bremsstrahlung of the ^''''Ce 
/?-decay electrons slowing down in the Cerium material. 



Compton scattering and photoelectric effect. The inten- 
sity of 2185 keV 7 rays is decreased by a factor < 10^^^ 
{Xw ~1.2 cm) [14|, to reach a tolerable rate. 

The energy spectrum of external bremsstrahlung pho- 
tons in the cerium is estimated with a simulation using 
the cross section of 13]. Results were confirmed with a 
GEANT4 [22] simulation. The number of photons above 
a prompt signal threshold of 0.9 MeV is 6.5 • 10^'^ photons 
per /? decay, and lO"** photon per /3 decay >2.0 MeV. 

An important remaining background source could 
be the W shield itself. Activities at the level of ten to 
hundreds mBq/kg have been reported. We anticipate 
the need of an external layer of ultrapure copper, set 
to 2 cm in our simulation. It allows one to achieve the 
radiopurity and material compatibility requirements. 
Assuming a ^4 tons shield we consider a prompt and 



delayed event rates of 50 Hz and 25 Hz, respectively. The 
escaping 7 are attenuated in the LS, assuming a 20 cm 
attenuation length [lj|. Beyond a distance of 1.5 m from 
the source the backgrounds become negligible. Any of 
the bremsstrahlung photons or shielding backgrounds 
can account for either the prompt or delayed event, 
depending on their energy. The sum of the backgrounds 
integrated over their energy spectrum is shown on 
Fig. [21 supporting the case of kCi Vf. source versus MCi 
j/g source for which solar-z^'s become an irreducible 
background. A light doping of the LS with gadolinium 
or an oil buffer surrounding the shielding would further 
suppress backgrounds; finally, non-source backgrounds 
could be measured in-situ during a blank run with an 
empty shielding. 

We now assess the sensitivity of an experiment with a 
50 kCi ^''^Ce source running for 1 year. With the shield 
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FIG. 4. 95% CL. exclusion limit of the 50 kCi.y "*Ce exper- 
iment obtained in the Am^ew and sin'^(2^„ow) plane (2 dof). 
Our result (red and green lines with and w/o knowledge of 
source activity) is compared to the 90% and 95% CL. inclu- 
sion domains given by the combination of reactor neutrino ex- 
periments, Gallex and Sage calibration sources experiments, 
MiniBooNE as described in Fig. 8 of _3] (blue lines). 



described above and using events between 1.5 m and 6 m, 
the background is negligible. With t^vr?^^^ = 2 eV^ and 
sin^(20„ew) = 0.1, the interaction rate decreases from 
40,000 to 38,000 per year. The 95% CL. sensitivity is 
extracted through the following function: 






X =2.2^ 



K^s 



a 



uJiV^^p 






(1 + a^A^^i^p) 



— , (3) 



where N^\^^ are the simulated data in the no-oscillation 
case and A"exp the expectations for a given oscillation 
scenario, in each energy Ei and radius Rj bin. ai, is 
a 2% fully uncorrelated systematic error, accounting 



for a fiducial volume uncertainty of 1% in a calibrated 
detector, as well as for (e"*", n) space-time coincidence 
detection efficiencies uncertainties at the sub-percent 
level. (Tjv is a normalization error of 1%, describing 
for the source activity uncertainty (from calorimetric 
measurement, see [21J), and a is the associated nuisance 
parameter. Fig. |4] clearly shows that 50 kCi of ^'^''Ce 
allows us to probe most of the reactor antineutrino 
anomaly parameter space [3| at 95% CL. An analysis 
assuming no knowledge on the source activity shows 
that the oscillatory behavior can be established for 
^"^ncw < 10 eV^. We note that a 10 kCi source would 
be enough to test the anomaly at 90% CL. 

The reactor D,, anomaly and the hypothetical existence 
of a 4*'' ly state mixing with the fe could be efficiently 
tested by deploying a /3-source of about 10 grams of ^^"^Ce 
(or ^°^Ru) at the center of a large low background liquid 
scintillator detector. The technical challenge lies in the 
production of the source itself, as well as in the realization 
of a thick ultra pure W/Cu shielding surrounding the Ce 
(Ru) material. Significant collaborative work would be 
needed to bring this idea to fruition. We thank A. lanni, 
M. Pallavicini, B. Littlejohn, and J. Rich for discussions, 
and support by the "Origin and Structure of the Uni- 
verse: Cluster of Excellence for Fundamental Physics" 
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